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Abstract 

If the gravitino G is very hght and ah the other super symmetric particles are above 
threshold, supersymmetry may still be found at colliders, by looking at processes 
with only gravitinos and ordinary particles in the final state. We compute here the 
cross-sections for some distinctive signals at hadron colliders: photon plus missing 
energy, induced by qq — > GG'y, and jet plus missing energy, induced by qq — > GGg, 
qg — > GGq, qg — > GGq and gg GGg. From the present Tevatron data, we estimate 
the bound > 2.3 x 10~^ eV on the gravitino mass, corresponding to the bound 
Vf > 310 GeV on the supersymmetry-breaking scale. We foresee that the upgraded 
Tevatron and the LHC will be sensitive to values of up to 4.0 x 10~^ eV and 
6.2 X 10~^ eV, corresponding to ^/F up to 410 GeV and 1.6 TeV, respectively. 
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contains important changes. The correction of a sign error modifies the relevant partonic 
cross- sections. The sensitivity to \fF (m^/2) is only slightly weakened. 
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1 Introduction 



While there are strong theoretical motivations for supersymmetry, with supersymmetry- 
breaking mass splittings of the order of the electroweak scale (for reviews and references, 
see e.g. there is at present no compelling argument to select a definite value of the 
supersymmetry-breaking scale ^JF or, equivalentlyQ, the gravitino mass 7/13/2. As reviewed 
in Pi , different possibilities should then be kept in mind when performing phenomenolog- 
ical analyses: they are characterized by a 'heavy', 'light' or 'superlight' gravitino. In the 
heavy gravitino case, reactions involving the gravitino are never important for collider 
physics. In the light gravitino case, the gravitino can be relevant in the decays of other su- 
persymmetric particles, if there is sufficient energy to produce the latter. In the superlight 
gravitino case, also the direct production of gravitinos (with or without other supersym- 
metric particles) can become relevant. In this paper, we concentrate on the superlight 
gravitino case, where ^J~F can be close to the electroweak scale and, correspondingly, 7/13/2 
can be several orders of magnitude below the eV scale. 

Many aspects of the superlight gravitino phenomenology at colliders have been dis- 
cussed long ago 0, and also more recently [Q. In all these papers, however, it was 
assumed that some other supersymmetric particle, for example a neutralino or one of the 
spin-0 partners of the gravitino, is light enough to be produced on-shell in some reaction. 
Here we take an orthogonal point of view: there may be experiments where the available 
energy is still insufficient for the on-shell production of other supersymmetric particles, 
but nevertheless sufficient to give rise to final states with only gravitinos and ordinary 
particles, at measurable rates. In a recent paper [Q, some of us considered the process 
e^e~ —>■ GG'-f, which may give rise to a distinctive photon + missing energy signal at e'^e" 
colliders. In the present work, we study the possible signals of a superlight gravitino at 
hadron colliders, such as the Tevatron or the LHC At the partonic level, we consider the 
following subprocesses: 



q + q + G + (1) 

q + q G + G + g, (2) 

q + g G + G + g, (3) 

q + g G + G + q, (4) 

g + g ^ G + G + g. (5) 



The process of eq. ([l|) corresponds to a photon + missing transverse energy (^t) sig- 
nal, those of eqs. (|D-(^ to a jet + signal. We compute the cross-section and the 
relevant angular distributions for the processes of eqs. (|l])-(^), in the limit in which the 
supersymmetric particles of the Minimal Supersymmetric Standard Model (MSSM), such 
as squarks and gauginos, and all other exotic particles, such as the spin-0 partners of the 
goldstino, are heavy. We then analyse the resulting phenomenology at the Tevatron (pp, 
^/S = 1.8 TeV, L ~ 100 pb"^), the upgraded Tevatron {pp, = 2 TeV, L ~ 2 fb"^), and 

^We recall that, in a flat space-time, F = \/3m3/2Afp, where Mp = {SnGNy^/^ ~ 2.4 x 10^® GeV is 
the Planck mass. Note that we take F real and positive, which is not restrictive for the present paper. 
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the LHC {pp, -v/S* = 14 TeV, L ~ 10 fb~^). We show how the present Tevatron data, in 
the absence of a signal over the Standard Model background, should allow us to establish 
the lower bound v'F > 310 GeV, or, equivalently, ^3/2 ^ 2.3 x 10~^ eV. These bounds are 
considerably better than the one estimated from the present LEP data, a/F > 200 GeV 
or ms/2 > 10~^ eV. We also estimate the future sensitivity of the upgraded Tevatron and 
the LHC: ^/F up to 410 GeV and 1.6 TeV, corresponding to up to 4.0 x 10~^ eV 
and 6.2 x 10""^ eV, respectively. In contrast with other collider bounds discussed in the 
literature , those discussed in the present paper cannot be evaded by modifying the mass 
spectrum of the other supersymmetric particles: making some additional supersymmetric 
particle light leads in general to stronger bounds. Therefore, making use of our results, 
an absolute lower bound on the gravitino mass can be established. 

The paper is organized as follows. In sect. 2, we spell out and discuss the low-energy 
effective Lagrangian used for our calculations. In sect. 3, we compute the relevant differ- 
ential cross-sections at the partonic level. In sects. 4 and 5, we discuss the phenomenology 
of the jet + and 7+ signals at hadron colliders, emphasizing the implications of the 
present Tevatron data and the prospects for the upgraded Tevatron and the LHC. Finally, 
in sect. 6 we summarize and discuss our results. 

2 The effective Lagrangian 

For the theoretically oriented readers, we summarize the framework of our calculation. 
We are considering processes whose typical energies are much larger than the gravitino 
mass, but smaller than the masses of all the other particles not belonging to the Standard 
Model. We can then work with the low-energy effective theory defined by the following 
two-step procedure. We start from a generic supergravity Lagrangian, assuming that su- 
persymmetry is spontaneously broken with a light gravitino, and we take the appropriate 
low-energy limit: Mp — > 00 with y/T" fixed. In accordance with the supersymmetric equiv- 
alence theorem 0, gravitational interactions are consistently neglected in this limit, and 
we end up with an effective (non-renormalizable) theory with linearly realized, although 
spontaneously broken, global supersymmetry, whose building blocks are the chiral and 
vector supermultiplets containing the light degrees of freedom, including the goldstino. To 
simplify the discussion, we assume that such a theory has pure F-breaking and negligible 
higher-derivative terms. Since we are interested in the case where the available energy is 
smaller than the supersymmetry-breaking mass splittings, we perform a second step and 
move to a 'more effective' theory, by explicitly integrating out the heavy superpartners 
in the low-energy limit. The only degrees of freedom left are then the goldstino and the 
Standard Model particles, and supersymmetry is non- linearly realized 0, H, . 

The derivation of the low-energy effective Lagrangian, to be used as the starting point 
for our calculations, proceeds exactly as in P|. The only differences are that now the 
unbroken gauge group is SU{3)c x U{l)em and that we are dealing with quarks instead of 
leptons. In analogy with we neglect the light quark masses and any possible mixing 
in the squark sector. Denoting with the light quark fields (where a = 1,2,3 and 
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q = u,d, s, c, b), and with (A^)^ = A^{T^)'^ the gluon fields (where A = 1, . . . , 8 and 
are the SU{3)c generators), the covariant derivative with respect to SU{3)c x U{l)em 
reads[|: 

(D^); = (d, + leQ.A^) 5^ + igsiA.Tp ■ (6) 
Focusing only on the terms relevant for our calculations, we obtain: 

^e// = -\f,.f^^ - -F^^F^^^ + iY.Umiq' + Ig^g + j:o,. (7) 

q i=l 

The terms Oi {i = 1, . . . ,4) are local operators, bilinear in the goldstino and generated by 
the exchange of massive particles in the large- mass limit. Their field- dependent parts have 
mass dimension d = 8, and the corresponding coefficients scale as l/F"^. The operator Oi 
involves only vector bosons (photons or gluons) and goldstinos, and is generated by 
integrating over the spin-0 partners of the goldstino and the gauginos: 

Oi = (g [7^ 7l F,^ ^ [Y, r] F,^G + G [7^ Y] F^AP)^^ [Y, Y] Ff„ g) , (8) 

where p is constructed with the covariant derivative in the adjoint representation. The 
operator O2 is a four-fermion interaction among quarks and goldstinos, recovered by com- 
bining a contact term in the original Lagrangian with contributions originating from squark 
exchanges 0: 

^2 = {^aGD (Cg") - g„75Gn (^75?") } ■ (9) 

The operator is generated by attaching a photon or a gluon to a squark exchanged 
among quarks and goldstinos^: 

^3 = E [eQ.^'^^l + 9s{A,r^ [q^Gd, (Cq^) - q^j.Gd, (^75?^) } • (10) 

The operators O2 and O3 are not gauge invariant, but their sum is contained in a gauge- 
invariant combination, as can easily be verified. Finally, the operator O4 is a gauge- 
invariant contact term that directly contributes (as O3) to the scattering amplitudes of 
eqs. (p-(^). It originates from a combined squark and gaugino exchange: 



8F2 



eQ,F,J^^ + gsF^^iT^)'^] {g [Y, Y] {q^q^G - l,q^q^l,G 



+ {GqPq^ - G^.q^q^l,) [7^ Y] g} . (11) 
As a final remark, we should mention that the processes considered in the present pa- 



per were also considered in |]T0[, in the same kinematical limit, but relying on the effective 
Lagrangian traditionally associated with the standard non-linear realization of supersym- 
metry 0, which was believed at that time to be unique at leading order. As was recently 
emphasized [Q, however, such an approach is not equivalent to ours, and the two lead 

^Notice the sign change with respect to the first version of this paper, following the analogous sign 
change in ref. 
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to different results, both consistent with supersymmetry. For a model-independent study, 
we would need the general form of the low-energy effective interactions, allowed by the 
non-linearly realized supersymmetry, that may contribute to the relevant amplitudes at 
leading order. This was not known until recently, when a definite theoretical prescription 
for such an investigation became available (the same paper also confirmed the results 
of [Q, and found the explicit general form of the coupling of two on-shell goldstinos to a 
single photon). We are postponing a general phenomenological analysis to a forthcoming 
paper. Because of the strong and universal power-law behaviour of the cross-sections, 
always proportional to s^/F^, we expect our results to be rather stable with respect to 
variations of the parameters characterizing the most general non-linear realization. How- 
ever, should a signal show up at the Tevatron or the LHC, having the general expression 
of the cross-section would be very important, since a detailed analysis of the spectrum and 
of the angular distributions of photons, jets and missing energy would offer the unique 
opportunity of distinguishing among possible fundamental theories. 



3 Partonic cross-sections 

We now compute the cross-sections for the processes of eqs. (|lD-(|D, and comment on the 
cross-section for the process of eq. (H). 

For the process qq —>■ GG'j, the calculation is a trivial extension of the one performed 
in for e"'"e~ GG'-f. We give here, for future use, a compact expression for the matrix 
element squared, summed over the helicities of the initial and final states. The calculation 
has been performed both by standard trace techniques, with the help of the program 
Tracer and by helicity-amplitude techniques |T2|. We denote by {pi,P2,qi,q2,k) the 



four-momenta of the incoming quark and antiquark and of the outgoing gravitinos and 
photon, respectively. For any given quark flavour, the result is 



(|M(gg^GG7)P) 



Q^a 327r 



X [^(pi,P2,gi,g2,fc) + A(p2,pi,gi,g2,fc) + (gi ^ 92)] , 

where the symbol (. . .) denotes the average (sum) over colours and helicities of the initial 
(final) state, Qq is the electric charge of the quark under consideration, a = e^/ (An) is the 
electromagnetic fine-structure constant, iV = 3 is the number of colours and 

A{pi,P2,qi,q2,k) = 

{kpi){kp2f{kqi){pip2){piq2) - {kpi){kp2){kqif{pip2){piq2) 
-{kpi){kp2f{kqi){piqi){piq2) + {kpi){kp2){kqif{piqi){piq2) 
-2{kpi){kp2){kqi){piP2){piq2f + {kpi){kqif{piP2){piq2f 
+ {kpi){kp2){kqi){piqi){piq2f + {kpif{kp2){kqi){piq2){p2qi) 
-{kp2){kq2){piP2){piqi){piq2){p2qi) + {kp2){kq2){piqif{piq2){p2qi) 
-{kpif{kqi){piq2f{p2qi) + {kp2){kqi){piP2){piq2)^{p2qi) 
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-{kp2){kqi){piqi){piq2f{P2qi) - {kpi){kq2){piqi){piq2){p2qif 

+ {kp2){kq2){piP2){P2qif - {kq2){PiP2){Piq2){P2qif 

+ {kpi){kqi){piqi){piq2){p2qi){p2q2) + {kqi){piP2){piq2){P2qif {P2 q2 ) 

+ {kpi){kp2){piP2){piq2){P2qi){qiq2) - {kpi){kqi){piP2){piq2){p2qi){qiq2) 

-{kpi){kp2){piqi){piq2){p2qi){qiq2) + {kpif{piq2){p2qi)^{qiq2) 

-{kp2){piP2){piq2){p2qif{qiq2) + {piP2){piq2f{p2qif{qiq2) ■ (13) 

After integrating over part of the phase space, we obtain: 



dx^d cosO^ N 1607r2F4 



■ /(x^,cos6'^) , (14) 



where (a/s, x^, 6^) are the total energy, the fraction of beam energy carried by the photon, 
and the photon scattering angle with respect to the quark, in the centre-of-mass frame of 
the incoming partons. The function / is given by [0 



/(x^, cos6'-y) = {1 — X. 



'l — x^){2 — 2x^ + x'^^) x^(— 6 + + x^) sin^ 6'^ 



+ 



x^ sin 9^ 16 32 



The extension to the process qq GGg is straightforward, since the only difference is 
the replacement of the photon with the gluon. By replacing a with Cpas in eqs. (|12D 
and dll), where Cp = {N"^ - l)/2N = 4/3, we obtain: 

da Cpas s . f^Xg^cosOg) . (16) 



dxgd cos Og N 1607r2F4 

For the crossed channels, qg — > GGq {qg — > GGq), analogous results can easily be 
obtained by taking eqs. (|12D and (|13]), by replacing Q'^a/N with Gpas/iN"^ — 1), by 
performing the substitution k —p2 {k ^ —Pi), and by changing the overall sign. 
After integrating over part of the phase space, we easily obtain the one-parton inclusive 
distributions in terms of the partonic centre-of-mass beam-energy fraction Xq{xq) and 
scattering angle 9g{9g). For the qg — > GGq channel: 

Gpas s gl^Xq^cosOq) , (17) 



dxgdcosOg 1 1607r2F4 



where 6q is the angle between the direction of the incoming and outgoing quarks, and 

l + xl (l-x,)(l-2x, + 2x2) 



g{Xq, COS 6^5) = (1 — X 



\2 



qJ 



16sin2(^^/2) 4cos2(^g/2) 

(1^ 



^ Xg(12-27xg + 12x^) ^ x^(3-4xg)cose, 



32 32 

Because of the obvious symmetry properties of eq. (|l^), an identical expression is obtained 
for the qg —>■ GGq channel. 
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PM , (21) 



It is easy to verify that all the above results satisfy the required factorization properties 
in the case of coUinear emission: 

^ GG,) * ^ GG; (1 - ;^J-L- P„(,,, , (19, 
GG,') ^4' GG; (1 - ^^j-^ P„(.,) . (20) 

. „ (19 ^ GGq') '-i" a(qq ^ GG; {1 - x,)s) ^ ^] 

dXqdcosOq 71 4cos^(6'g/2) 

where are the standard Altarelli-Parisi splitting kernels, and 

1 

a(qq^GG:s) = (22) 
- GG; .) = ^ ^ . (23) 

We did not perform the complete calculation of the gg GGg process, since the con- 
tribution of the gg channel is significantly suppressed at the highest values of s, where 
the signal is more hkely to emerge over the QCD backgrounds. To verify this statement, 
we shall use in the next section an approximate estimate of the gg process, based on the 
coUinear approximation: 

_^(,, ^ GG,) ~ a(,, ^ GG; (1 - ^ -L_ P,,(.J . (24) 

Our results will confirm that indeed the contribution of this channel is strongly suppressed 
in the kinematical regions of interest. 



4 The jet-plus-missing-energy signal 

In this section, we evaluate the rates for the production of a single jet of large transverse 
energy {Et) plus a gravitino pair in hadronic collisions. We shall consider the cases of 
the currently available Tevatron data (TeV I, pp collisions at \fS — 1.8 TeV, with an 
integrated luminosity of approximately 100 pb^^), of the upcoming upgraded Tevatron 
run (TeV II, pp collisions at \fS = 2 TeV, with an integrated luminosity of approximately 
2 fb~^), and of the LHC {pp collisions at = 14 TeV, with an integrated luminosity 
of approximately lOfb^^). The dominant irreducible background to the signal comes 
from the associated production of a Z vector boson and a jet, with the Z decaying to a 
neutrino pair. Given the large values of Et we will be considering (in excess of 100 GeV), 
no significant detector backgrounds are expected. Notice that the standard experimental 
searches of supersymmetry using multi-jet-plus-^T final states require the presence of at 
least three jets. No experimental study has been published, to date, of the single-jet-plus- 
f^T signal discussed here. As a result we have to rely on our own background estimate, 
since nothing is available from the experimental studies. We shall therefore extract the 
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Figure 1: Production rates for pp GG+jet and pp jet+i^z/ at VS = 1.8 TeV. 



F from the comparison of the expected 



region of sensitivity to the value of the scale 
signal and the calculated background rates. 

In our numerical analysis, we shall use the parton-density parametrization MRSR2 [|13 
with a renormalization/factorization scale fi = Et- Since the cross-sections for both signal 
and background processes scale linearly with as, their ratio is rather insensitive to this 
choice. In particular, since the signal rate is proportional to F~^, the limits we shall set on 
\/F are not significantly affected by this choice, as we have explicitly verified by varying 
the scale /i in the range Et/4: < /i < AEt, independently for signal and background. 

Figures show the contribution to the signal rate of the three possible initial states 
qq, qg + gq and gg. The distributions shown correspond to the jet rate integrated over 
transverse energies larger than a given threshold E^^"', and within a jet-rapidity range 
I77I < 2.5. For the three figures we selected the values of \/F = 310, 410, and 1600 GeV, 
respectively. This particular choice will be justified later on. The relative contributions of 
the three channels and the shapes of their distributions are not affected by the value of \/F. 

For comparison, we also show the background rate for the process p'p^ Z jet^ uu jet. 
This was calculated at leading order in QCD. The capability of LO QCD to describe this 
process is proved by the recent study of the Z+ multijet final states published by CDF |]T| 



with the Z decaying to charged lepton pairs. In this paper the ratio Data/Theory for the 
Z-plus-l-jet final state was found to be 1.29 ± 0.17, using /i = Et- For consistency with 
this result, we rescaled our LO background prediction by a factor 1.29. To be conservative, 
however, we did not rescale the signal rates. 

We emphasize the following features of these distributions. First of all, notice that 
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Figure 2: Production rates for pp — > GG + jet and pp — > jet + vu at VS = 2 TeV. 
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Figure 3: Production rates for pp — > GG+jet and pp jet + z/P at y/S = U TeV. 
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Figure 4: Solid line (scale on the left axis): estimated 95% CL exclusion limits 
on \/F from the absence of a deviation from the SM jet+^r rate in events with 
jet-£'r > EJp^^ at TeV I. Dashed line (scale on the right axis): number of jet + GG 
events with jet-i?T > -E™" corresponding to the 95% CL value of ^/F. 

the shape of the background is very similar to that of the signal. This suggests that the 
best limit on \/F will be extracted from a comparison of signal and background in the 
widest possible region of Et compatible with other experimental requirements. Secondly, 
notice that the contribution of the gg initial state is always negligible in the interesting 
regions of Et-, due to the smallness of the gluon density compared with the valence quark 
density. This is true regardless of the colhder energy and justifies our use of an approximate 
expression for the gg GGg matrix elements. Finally, notice that while at the Tevatron 
the signal is dominated by the qq initial state, the dominant contribution at the LHC 
comes from the qg + gq channel. This is because the antiquark density in pp collisions is 
much smaller than the gluon density. 

To extract the region of sensitivity to the scale a/F, we count the number of background 
events expected above a given threshold ii^™*" and evaluate the value of \fF that can 
be excluded at the 95% confidence level (CL). For values of i?™" such that no more 
background events are expected, this corresponds to a limit of 3 signal events. The 95% CL 
lower limits on \/F for the three cases considered are plotted as functions of i?™" in 
figs. ^-|^. In the same figures we also plotted (as dashed lines) the number of signal events 
corresponding to the selected value of \fF. The best limits on a/F that can be extracted 
with the chosen integrated luminosities are given by 310, 410 and 1600 GeV for TeV I, 
TeV II and the LHC, respectively. The corresponding limits on are 2.3 x 10~^ eV, 
4.0 X 10"^ eV and 6.2 x 10"^ eV. These limits are obtained for E^*" = 100 GeV in the 
case of TeV I and TeV II, and for EJp"^^ = 200 GeV in the case of the LHC. Lower values 
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Figure 7: Production rate for pp GG + 7 at \fS = 1.8 TeV. 



of i?™" could lead to more stringent limits, but in this case other backgrounds could 
become important, and a detailed experimental analysis would be required. We should 
also mention that our procedure relies on the detailed prediction of the background rates. 
Notice nevertheless that in the regions of Et where a large background is expected, its 
absolute normalization can be determined with high accuracy by measuring the [Z 

+ jet rates. In these regions of large rates, furthermore, one could use not just the 
total rate above a given threshold, but the ii^T-distribution as well. 



5 The photon-plus-missing-energy signal 

While the production rate of events with a photon plus the gravitino pair is smaller than 
that with a jet, the relevant source of background (associated production of ■Z'7) is also 
smaller. A recent analysis from DO |]l5l, for example, did not find any 7+ event with 
Ex'y > 70 GeV in approximately 13 pb"^ of data. It is therefore interesting to examine 
the constraints on ^/F set by searches in this channel. Figure ^ shows the 7+ signal 
rate, as a function of the photon Et threshold. We chose for this plot a/F = 210 GeV. In 
the case of photon final states, only the qq annihilation channel contributes. The results 
we find agree with the naive estimate: 

a{GG^) ^^^3^^_ (25) 
ct(G'G jet) Cpas 

The 95% CL exclusion limit on ^/F as a function of i?™", for the 13 pb^^ of data cor- 
responding to the available DO analysis, is given in fig. ^. In addition to imposing a 
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Figure 8: 95% CL exclusion limits on \'F from the absence of 7+ events at 
DO 0. 

geometrical acceptance cut of \ri^\ < 2.5, we reduced both signal and background rates by 
the photon identification efficiency = 0.6 [ITsj. For E'™" = 70 GeV, the threshold value 
above which DO reports no events found, we obtain ^/F > 210 GeV. Slightly better con- 
straints could be obtained by using lower thresholds. We look forward to the completion 
of the analyses of the full 100 pb~^ current data sample, for which we estimate a 95% CL 
limit of 250 GeV (265 GeV) for E^'" = 70 GeV (50 GeV). These numbers are only a 
factor of 15% worse than those obtained with the jet+^-r analyses. Further improvements 
may arise from a combination of the DO and CDF statistics. 

The sensitivity achievable in the 7+ channel at TeV II is displayed in figs. |^ and 
rUl which show the production rates and 95% CL limits after accumulation of 2 fb~^ of 



data at 2 TeV. Similar plots for the LHC are shown in figs. |TT] and |T^. In this case we set 
the photon-identification efficiency equal to 1, in the absence of a determination based on 
real data. While the Tevatron sensitivity in the photon channel is almost comparable to 
that in the jet channel, a clear advantage of the jet signal over the photon one is seen at 
the LHC. This is due to the smaller luminosity of the qq initial state relative to the qg 
initial state at the LHC. 



6 Discussion and conclusions 

We presented in this paper a phenomenological analysis of high-energy collider constraints 
on the scale of supersymmetry breaking ^/F, in one-to-one correspondence with the grav- 
itino mass ^3/2, in models where the only light supersymmetric particle is the gravitino. 
A typical signature for these models is the annihilation of e'^e~ or qq pairs into gravitino 
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Figure 10: Solid line (scale on the left axis): estimated 95% CL exclusion limits 
on s/F from the absence of a deviation from the SM 7+ fix rate in events with 
^-Et > -Ej?'*" at TeV II. Dashed line (scale on the right axis): number of GG + 7 
events with j-Et > -Ej?'*" corresponding to the 95% CL value of s/F. 
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Figure 12: Solid line (scale on the left axis): estimated 95% CL exclusion limits 
on s/F from the absence of a deviation from the SM 7+ rate in events with 
j-Et > -E'j?'*" at the LHC. Dashed line (scale on the right axis): number of GG+j 
events with ^-Et > -Ej?*" corresponding to the 95% CL value of VF. 
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Figure 13: y/s distribution for events with Et > 100 GeV (solid line) and 
Et > 200 GeV (dashed line) at the Tevatron. 



pairs. To tag the events, additional emission of either a jet or a photon is necessary. 
The study of LEP data in ref. has been extended here to signals from the Tevatron 
collider. The limits from the Tevatron are significantly stronger than those obtained 
from LEP, thanks to the higher energy reach. Supersymmetry-breaking scales as large as 
\/F = 210 GeV (m3/2 = 1.1 x 10~^ eV) are already excluded at 95% CL by a DO search 
for events with photons and missing energy. The analyses of the full TeV I data sample 
should increase this limit to ^/F ~ 265 GeV (^3/2 — 1.7 x 10~^ eV). Comparable and 
probably better limits should be obtained from the current data in the jet plus final 
state, where we foresee sensitivity up to — 310 GeV (m3/2 — 2.3 x 10^^ eV). The 
upgraded Tevatron, accumulating 2 fb^^ at 2 TeV after the year 2000, should improve the 
sensitivity up to \/F ~ 410 GeV (m3/2 — 4.0 x 10~^ eV). The sensitivity will then increase 
to ^ 1.6 TeV (m3/2 — 6.2 x 10~^ eV) once the LHC data become available. 

These estimates have been obtained by using an effective Lagrangian containing only 
the degrees of freedom associated to the ordinary particles and the ±1/2 helicity states of 
the gravitino. Such a Lagrangian is non-renormalizable and is expected to provide a good 
approximation of the fundamental theory in an energy range bounded from above by a 
critical value E^, ~ (2-^3) v'F. Beyond Ec, perturbative unitarity gets violated and the non- 
local character of the low-energy theory should show up. These considerations also bound 
the mass M of the new particles called to rescue unitarity: M < (2 ^ 3)^/F. We should 
therefore check that, for each reaction considered, the typical energy ^/s probed by the 
partonic process remains within the allowed domain. Similarly, to be consistent, we should 
be able to consider superpartner masses sufficiently large to suppress direct production or 
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large threshold effects, while remaining within a few times the supersymmetry breaking 
scale /F g. 

These constraints turn out to hold for the cases considered in this paper. As an 
example, we show in fig. |1^ the distribution of the variable i/s for jet+^T events with 
j^mm _ ]^gg 200 GcV at the Tevatron. The figure shows that, for the typical sensitivity 
range of \fF ~ 300 GeV, the average value of ^/s probed by this class of events is indeed 
of the order of (2 ^ 3) \fF. We verified that similar results hold at the LHC, with average 
values of a/s in the range of 3 TeV. 

The study presented in this paper should allow experimentalists to establish absolute 
lower bounds on \/F and 7713/2- Indeed, we do not expect that our bounds dissolve as the 
threshold for superparticle production is approached or crossed. On the contrary, in that 
case possibly stronger limits on the gravitino mass could be put by exploiting the pro- 
duction of the gravitino in association with squarks and gluinos, or the production of jets 
in association with the spin-0 partners of the goldstino. No conclusive phenomenological 
study of this type has however been carried out as yet, since no experimental analyses are 
available of final states with one jet plus ^t- 

The reader may have noticed that, in our analysis, all weak- interaction effects were 
neglected. For the signals discussed, and in the case of a goldstino neutral under the full 
electroweak gauge group, this is a consistent approximation: virtual W and Z exchange 
cannot give appreciable contributions. We may ask, however, if the associated production 
of a gravitino pair and a weak vector boson may lead to detectable signals. For the GGZ 
final state, the Standard Model background [ZZ production with one of the Z decaying 
to a neutrino pair) is very small: only a few ZZ events have been observed up to now 
at TeV I. The signal, however, is also expected to be strongly suppressed with respect to 
the photon case, since the gain due to the coupling constant enhancement is much less 
than the loss due to the small leptonic branching ratio (hadronic decays suffer from severe 
background problems) and to the effects of the finite Z mass. For the GGW final state 
the background is large, being dominated by production of a single off-shell W . In both 
cases, we expect weaker limits on \/F and 7/13/2 than the ones discussed in this paper. 

We may also ask if the process 57 — > GGq, analogous to the ones considered in the 
present paper, may give rise to detectable signals at HERA. We studied this process in 
ep colhsions at ^/s = 300 GeV. With 100 pb~^, roughly corresponding to the present 
integrated luminosity collected by the ZEUS and HI experiments, we found no sensitivity 
at HERA for values of ^/F in excess of 100 GeV. 

In conclusion, we expect that hadron colliders will lead the search for a superlight 
gravitino in the next decade. 

■^We recall that the complete absence of supersymmetric particles with masses below the TeV scale 
is generically disfavoured by naturalness considerations, even if moderate exceptions may exist. For this 



reason, the limit on VF obtained from our study of the future LHC data essentially saturates the range of 
applicability of our formalism. Note also that our approach does not set interesting constraints on models 
with gauge-mediated supersymmetry breaking (for a review see ref. |16 ), in which typical superpartner 
masses are much smaller than VF, so that one is lead to set \/F > 10 — 100 TeV. 
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